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Dipole fluid as a basic model for the equation of state of ionic liquid 
in the vicinity of their critical point 
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The model of dipole fluid for the ionic liquids similar to the molten NaCl is proposed. The 
, estimates for the critical parameters are obtained with the help of the van der Waals equation of 

■ state. The influence of the rotation on the characteristics of a dipole pair and the location of the 

critical point is discussed. The dissociation of such fluid near the critical point is considered. 
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Introduction 

There is the strong indication for existence of the liquid-gas critical point in ionic liquids like molten NaCl. The main 
model for such systems is called Primitive Model or its restricted version where the characteristics (mass, diameter 
etc.) of the opposite charges are the same. Because of formation of associated ionic pairs the derivation of mean 
field equation of state in rigorous thermodynamic approach is yet unsolved. The existence of additional constraints, 
so called sum rules [jj], makes the problem of finding thermodynamically consistent mean field theory hard to solve 
comparing with the molecular liquid case. The various approaches have been proposed recently for the derivation of 
equation of state (EOS) for the RPM. Mostly they are based either on the improvements of classical Debye-Hiickel 
(DH) Q and Bjerrum || theory of association or Mean Spherical Approximation (MSA) H taking into account the 
association and various effective interactions. For the reviews of the results see || 0, g, |j. 

The location of the critical point for RPM varies in dependence of the interactions included into free energy. It is 
accepted that the following dimensionless parameters for the temperature T, particle number density n and pressure 
' P are used: 

<£; T* = ^, p* = na 3 , P* = ^j- A (1) 

CO ! 

, where a is the ionic diameter q is the absolute value of the charge. The characteristic value for T* takes the values in 
the interval 0.05 < T* < 0.08. The situation with the critical density p c is less definite. Its value varies from p c = 0.02 
to p c = 0.08. It should be noted that most of the analytical results predict low density (p c = 0.026) and "high 
temperature" (T* = 0.06 — 0.08) critical point. It cannot be excluded that such low value for p* is connected with 
small association taken into account with the help of thermodynamic perturbation theory. However the reliability of 
such type estimates is not quite clear. Therefore the development of the alternative approach grounded on the dipole 
fluid model seem to be expedient. 

The Monte Carlo (MC) simulation studies jl(| [ll], [l2| performed in recent years give "high density" p c = 0.07 — 0.08 
and low temperature T* = 0.04 — 0.05 critical point. 

In the present paper we estimate the critical characteristics of the ionic molten in the framework of the dipole 
O ■ fluid approach. We will also show that the change of molecular parameters due to rotations influences essentially the 
' location of the critical point. 



I. QUALITATIVE ANALYSIS OF THE CRITICAL BEHAVIOR OF THE DIPOLE LIQUID 

_ _ i 

NaCl is the simplest example of ionic liquid. In solid state it is ionic crystal. Above the melting point the positions of 
ions become unfixed, but this liquid remains strongly dissociated. Due to thermal expansion at increasing temperature 
the dissociation degree diminishes and the molten salt passes to a dipole liquid. At further increasing temperature 
and decreasing the molten density the dissociation degree grows again and the molten salt becomes completely ionized 
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system. Thus, at some temperature (T±,T2) and density (ni,n2) intervals the molten NaCl can be considered as a 
dipolc liquid. Further the additional arguments for this assumption will be given. 

Let us consider the general properties of this liquid and, in particular, its main critical parameters: the critical 
temperature and density. The interparticle interaction in dipole system takes the form: 

U(l,2) = U w (l,2) + U dd (l,2) (2) 

where the first term describes the Van-der-Waals interaction between molecules and 

U dd = -\- - 3 (d> H fe) ) (3) 

the proper dipole-dipole interaction, m = \r± — r2 1 is the interparticle spacing. Note that the dipole moment of a 
pair is some function of the equilibrium distance between ions in a pair: 

r (<) _ r W 

di = d(r d ) ni , n+ = w 

| r w _ r w| 

i = 1, 2, e is dielectric permittivity. 

Because of the dipole-dipole interactions are relatively weak the angular distribution of dipole moments di is close 
to isotropic one. More exactly we assume that the two particle distribution function g{d\,d,2) can be approximated 
by the first two terms in the expansion: 

g{d x ,d 2 ) = 1-0U M (1, 2) + ..., 0=rV ( 4 ) 

The approximation of such a kind allows us to exclude the orientational degrees of freedom in the configurational 
integral with the help of the perturbation theory. In fact this procedure is equivalent to the usage of the isotropic 
potential 



U(r 12 ) =< 17(1,2) >= U w {r 12 ) - U d (f^j 

TT , — 2 R «d 2 » 2 
Ud - oP 



(5) 



Here a sa a + + a_ sa 2a, a + and a_ are the diameters of ions Na and CI correspondingly and for simplicity we put 
a + = a_ and neglect the difference in masses of the ions, << ... >> denotes the average with internal partition 

function of a pair. It is easy to check that the inequality |J7 u; (r 1 2)| << U d (777) takes place at all r 12 . Therefore 

further the contribution \U W \ will be ignored. 

It is essential, that the averaging procedure restricts the applicability region of the potential (^) by interparticle 
spacings a < r 12 which gives the size of the "averaged" dipole of order 2a. Though the value of a may be slightly 
less than 2a since rotating dipoles are not the same as hard spheres of diameter 2a. It is quite clear in view of 
scattering cross section for the hard rotating dumbells. At this level a should be considered as the parameter (in 
general temperature dependent) of the dipole-dipole potential. The procedure of its fixing in the critical point will be 
discussed further. 

To describe the properties of the molten NaCl within the interval, where it can be considered as a dipole liquid, we 
can use the potential with hard wall: 

( OO , j-12 < cr 



Such a potential leads to the Van-der-Waals equation of state 

ndksT A frrn , 



where 



P=T^ b -MT)ni (7) 



. U I'M rr / U{r, T)r 2 dr = ^-U d , b = |V (8) 
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and 



and rid is the pair number density. Therefore the overall density is n = 2nd- In dimensionless form (0) and (||) read 
as: 

_ _?t*_ _ icn 2 

2 -bp* 4 P ' W 

2-7T 

^( T *) = ^3«^» 2 (io) 

Here all spatial parameters is given in units of a. The value of parameter << >> depends on the internal structure 
of the pair. Though by the order of magnitude «rd >>~ 1, nevertheless from []l3| it follows that the characteristic 
values of dipole moments correspond to << >> < 1. 

Since a is connected with the size of the pair we model its temperature dependence via relation: 

a =« r d » 5 (11) 

where 5 is the fitting parameter. 

Note that a is temperature dependent which we assume the same as that for << r >>. The functions << r" >> 
will be determined below. Note that the vapor phase contacting with the liquid one is the gas of dipole molecules. 
The van der Waals EOS is appropriate approximation for EOS for such vapor phase. Thus we can get the the critical 
parameters of this system using the van der Waals theory of the critical point. 

The equation (Q) leads to the following equations for the critical temperature and density (note that rid — ri/2, 
where n is the total number density): 

T *c = ^§«r 2 d », Pi = -ip (12) 

The estimates for these parameters are straightforward if we put er = 2, and take into account that due to small dipole 
moment of NaCl << r 2 d >>— 1 (in units of a): 

T* = f «0.04, pS = ^«0.04, ^ = ^3-10- (13) 



P* 3 
Z c = = — ra 0.19 

PIT* 16 

which are in satisfactory agreement with the values 

T* = 0.055 , p* = 0.026 P* = 3.6 • 10" 4 , (14) 



Z c = 0.25 

obtained within extended Debye-Hiickel-Bjerrum theory & augmented with ion-dipole interaction. Our value of the 
critical density is greater due to the neglecting the dissociation of the dipole pairs. 

Now we need to consider the dipole pair as itself since the parameters of the potential (^j actually are the averages 
over the internal partition function of a pair and therefore are the temperature dependent functions. 



II. THE DISSOCIATION OF THE ROTATING DIPOLE LIQUID 

In previous section the model of undissociated ionic liquid consisting of rotating dipoles has been introduced. Here 
we investigate the internal structure of the bound pair of ions. We take into account the fact that the energy of 
interaction of a pair should include centrifugal energy together with Coulombic potential as in standard problem of 
two bodies interacting via central field. 

The dissociation temperature for NaCl-molten is determined by the effective potential of an ion within a rotating 
dipole which includes the centrifugal energy: 

kT d w -U e „ = y - ^ (15) 
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where I = fir 2 is the moment of inertia of the charge with reduced mass fi = ™^™ m = m/2. At such high 
temperatures all degrees of freedom are all in equilibrium and we can use the estimate 

E rot =(^j\=kT. (16) 

Note that the equilibrium distance between ions in a pair, which is determined by the minimum of effective potential 
(H) with the help of (|l6|), is 

a eq = ^> (17) 

which is exactly the Bjerrum size of the pair i? Bj J3j (see also fit). The choice of dl7| ) as the size of the ionic pair is 
inappropriate from the physical point of view at low temperatures T* << 1 fij|. It is natural that with lowering T the 
size of a pair should become smaller tending to a at T — > 0. That is why it was suggested to use it for 1/T* > 2 only. 

Let us consider this question within the picture formulated above. To be more correct we will include the rotational 
energy into association constant, which is proportional to the internal partition function of the pair |(| g]: 

R 

K{T*,R)=4:iT J exp(-(3U eff )r 2 dr (18) 

a 

In 2D case one can put R — oo because of the logarithmic growth of the electrostatic potential and get the estimation 
of Kosterlitz-Thouless temperature of dissociation |u| . In 3D case there is the problem with upper cutoff in such an 
approach where the association constant is identified with internal partition function of the ionic pair. 

To define the size of a pair following Bjerrum we investigate the extremal points of the integrand in (|l8|). Doing so 
we get two solutions: 



where 



A = — — — , J = fia . 

Here R+ is solution of the Bjerrum type (minimum of the integrand in (|l8|)), which as has been said above is 
inappropriate. i?_ is another solution corresponding to the maximum of the integrand, which has quite reasonable 
values and correct behavior at low T* . It is easy to check that asymptotically for low values of the temperature 
T* the value of K{T*) is formed mainly by the maximum of the integrand. In addition the appropriate limiting 
behavior to the hard-core contact at formal limit T* — > is hold provided that A = 1/2. This value of A is in full 
accordance with the virial theorem jl5) . All these facts confirms that we can treat the quantity R_ as the size of the 
pair even at "high" temperatures T* < 0.1. In addition i?_ never exceeds 2, i.e. the interparticle the distance when 
the influence of other pairs and charges on the effective potential can be neglected (see Fig. [l]). So we put i?_ as the 
physical cutoff for ((ll). The existence of such length scale was noted in fij] basing on the the numerical analysis of 
the function K(T* ,K)/ K(T* , R B: >), though only Coulomb potential was included in Boltzmann factor. It gives the 
rate at which K(T*,R) rises very rapidly to its plateau value. In our case we find the same behavior of K(T*, R) at 
small temperatures, T* < 0.04 (see Fig. ||). Finally we see that there is the natural temperature interval for dipole 
fluid which is bounded from above by the temperature 

T-upper — o (20) 



Therefore for T < T* pper the dipole is stable as itself. Note that the existence of the temperature (|2fj) reminds 
Kosterlitz-Thouless transition in 2D case. In particular, the centrifugal energy introduced above plays the role 
analogous to the chemical potential ". . . required to create a pair of particles of equal and opposite charge at a 



distance ro apart ..." in Kosterlitz-Thouless model |14 . However, in contrast to 2D case in 3D case there is no any 
divergence in the size of the pair and therefore in its polarizability, but its derivative on the temperature has singular 
behavior. In other words the temperature derivative of the polarizability is singular but not the polarizability itself. 
This inference might seem as mere an artifact of introducing the upper cut off in (|l8|). But it should be noted that 
taking dielectric permittivity e as the order parameter, which is directly connected with polarizability, we get exactly 
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FIG. 1: The temperature dependence of R- 
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FIG. 2: The ratio K(T* , R+)/K(T* ,5 ■ R-) as a function of <5 and T*. 



the divergence of its temperature derivative even in mean field approximation. This may serve as additional support 
for the conjecture of intensive breaking of the dipole pairs at the critical point observed in some numerical experiments 

The existence of the interaction between dipoles and the free charges provides additional instability mechanism for 
their dissociation thus reducing the temperature of " ideal" dissociation (|2(]) because of the polarization of the dipole 
in the external field of dipole-dipole potential (J6j) and Coulombic field of free charges. The consideration given above 
states that there are two characteristic transition in the dipolc-dipole fluid: 1) "dipole liquid - dipole gas" critical 
point of van der Waals type; 2) the smeared dissociation "transition" from associated state to almost completely 
dissociated one. This smeared transition can be characterized by the temperature on the binodal at which the degree 
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of dissociation is 1/2. 

The interaction between translational degrees of freedom of the ions is characterized by the critical temperature 
of the liquid-gas transition, while internal, rotational degrees of freedom are involved into dissociation of such dipole 
fluid. These degrees of freedom will strongly interact if the corresponding potentials will be of the same magnitude 
i.e. T c w Td- The additional confirmation of coincidence of such transitions is the high degree of dissociation above 
the critical point observed in numerical experiments jlo[ and theoretical models || , which incorporate dielectric 
permittivity resulting from the existence of the dipole pairs. This means that T d cannot be less than T c . All said 
above means that in such situation we can not rely on (0) since appropriate EOS should incorporate all relevant 
interactions which lead to phase separation. In particular the critical temperature is sensitive to the temperature 
dependence of the parameter A(T). 

The condition for the dissociation of a pair in external electrostatic field is: 

< d- E >= - « U eff » , (21) 

where 

d = d + d md (22) 
is the dipole moment, which consists of proper and induced dipole moments. Obviously, 

< d ■ E >= , < d md • E >= x < E 2 > (23) 
Here x is the polarizability of a pair. The averaging over the thermal equilibrium gives: 

<E 2 >=^, (24) 

where k is the dielectric susceptibility of the medium and V p h is the characteristic volume . It is connected with the 
polarizability of the cavity. According to the definition: 

k = e —^- = -xp* (25) 

So we get the equation for the temperature in dimensionless form: 

T* = -«U eff » (26) 



P*V ph 



Within the proposed approach we put V p h = ^-1% where l c is the radius of first coordination sphere. This is the 
minimal volume for which the conception of continuity of the medium can be applied. By the order of magnitude 
l c w 1.5 a. The solution of ( f26| ) gives the dependence T*{5). In order to fix the value of 6 in the critical point which 
determine the size of the pair we should equate T*(S) and T*(S) obtained above. This way we get: 

T* = 0.048 , p* = 0.054 , P* = 4.8 • 10~ 4 , a = 1.8 , Z c = 0.19, (27) 

which are close to those obtained above (p"3|). In notations of ||, a = 2ci2- In this work the estimate for the parameter 
ci2 from simple geometric considerations was given: 0.825 < 02 < 1.565. Thus our estimate is in this interval. From 
the results obtained above we can infer that the dipole fluid of rotating dipoles in the vicinity of its liquid-gas critical 
point is about to dissociate. Sure our consideration is incomplete since it does not take into account the existence of 
free charges. 

Finally we estimate the Ginzburg number by the formula used for the molecular liquids |lq] : 

\ 6 



where tq —<< r >>« a is the interparticle spacing within ionic pair and £0 is the amplitude of the correlation length 
for density fluctuations. Since the density fluctuations are connected with the ones for dipole pairs we put it to be 
equal £ > a. Using the parameters of the critical point found in (^7|) we get the estimate: 

Gi < 0.04 (29) 
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Discussion 



The ionic and dipolc liquids form two natural approximations to describe the critical properties of the systems 
similar the molten NaCl. In our paper we have estimated the main critical parameters for liquid with hard dipole as 
well as consider the influence of the effects arising due to softness of a dipole molecule. In particular the last is very 
important to describe the dielectric properties of a system near the critical point. Besides, the variation of molecule 
parameters due to the rotations allows us to determine the equilibrium size of a ionic pair. 

It is not excluded that the quantum corrections to internal states of the dipole pairs will also slightly change the 
estimates. In particular the temperature dependence of the vibrational contributions to the heat capacity can also 
be studied. The following step is to construct the equation of state for small " soft" dipole molecules and to take into 
account the dissociation process with the help of perturbation theory. The combination of such an approach with that 
developed in |17[ on the basis of ionic liquid allows to narrow the region of the most probable values for the critical 
parameters. 

Our estimate for the critical temperature correlates with the known analytical results. Note that most of the 
analytical approaches based on EOS for low density Coulombic system (DH, MSA etc.) where the dissociation is 
taken into account perturbatively. 

Within the dipole liquid approach we have obtained the estimate for the Ginzburg temperature and have shown 
that it less than one for the simple liquid by a factor 10 -2 — 10 _1 . The approximation of the dipole liquid allows us 
to analyze in the evident form the contribution of the polarizational effects |l7j . One can show that the lasts lead to 
the further considerable decrease of the Ginzburg temperature. 

Note also the possibility for the appearance of new inhomogeneous phase near the critical point of ionic liquids. 
Since the dissociation temperature Td is near T c , the system can desintegrate on the regions with the essentially 
different values of the ionization degree A: the drops of ionic and dipole liquids. As a consequence the region with the 
Ising-like behavior cannot be reached. This scenario needs in very careful investigation. These and other questions 
will be the subjects of further works. 
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